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WAVEGUIDE-MODE  POWER  BUDGET  FOR  AN 
ELF/VLF  TRANSMITTING  SATELLITE 


INTRODUCTION 

The  study  of  electromagnetic  propagation  fr^m  a satellite  to  ground  at  extremely 
low  frequencies  (ELF)  or  very  low  frequencies  (VLF)  is  an  interesting  and  important 
problem  since  LF  electromagnetic  waves  can  penetrate  to  significant  depths  below  the 
earth’s  surface  for  reception  at  submerged  or  buried  antennas.  These  waves  have  unusual 
characteristics  because  the  ionospheric  plasma  and  the  earth’s  magnetic  field  govern  the 
directions  in  which  the  waves  are  generated  and  propagated.  The  waves  exhibit  a strong 
tendency  to  fodow  the  direction  of  the  earth’s  magnetic  field  while  they  are  within  the 
ionosphere  (Fig.  1).  When  they  emerge  below  the  ionosphere,  they  enter  a leaky  wave- 
guide formed  between  the  edge  of  the  lower  ionosphere  and  the  earth’s  surface,  and 
propagate  via  waveguide  modes  to  be  received  at  distant  points  on  the  earth’s  surface.  In 
the  literature  this  waveguide  is  often  referred  to  as  the  earth-ionosphere  waveguide. 

Numerous  calculations  and  experiments  have  been  performed  concerning  the  propa- 
gation of  ELF  and  VLF  electromagnetic  waves  from  point  to  point  within  the  waveguide. 
Galejs  [1,2],  Einaudi  and  Wait  [3,4]  and  Pappert  [5]  have  previously  treated  the 
problem  of  the  excitation  of  waveguide  modes  by  transmitting  electric  or  magnetic 
dipoles  at  satellite  heights.  While  these  treatments  specified  the  magnetic  or  electric  dipole 
moment  required  at  the  satellite  to  achieve  a given  field  at  some  distant  point  on  the 
earth’s  surface,  they  did  not  further  specify  the  electrical  powers  that  are  required  to  pro- 
duce the  electric  or  ma  ;netic  dipole  moment  in  a given  realizable  antenna.  The  present 
report  adds  this  important  information  to  the  problem  for  the  transmitting  magnetic 
dipole  antenna.  Calculations  of  radiation  resistance  are  performed  for  various  transmitting 
antenna  sizes  and  heights  at  several  frequencies  using  the  methods  described  by  Baker 
[6] . The  ohmic  resistance  of  a loop  is  readily  calculable  from  the  loop  antenna  size 
parameters  knowing  the  conductivity  of  the  loop  material.  In  this  report  the  loop  antenna 
is  considered  to  be  made  of  aluminum  and  to  operate  at  room  temperature. 

The  transmitter  power  required  to  produce  a given  signal-to-noise  ratio  (S/N)  at  a 
given  point  on  the  earth’s  surface  is  strongly  dependent  on  the  exact  ciicomstances  under 
which  the  communications  circuit  is  needed.  For  example,  the  required  power  depends 
on  the  satellite  antenna  size,  altitude  and  composition,  the  transmission  frequency,  the 
required  S/N  ratio  for  communications,  the  time  of  day,  season  of  the  year  and  the  size 
and  geographical  location  of  the  desired  coverage  area.  The  first  results  presented  in  this 
paper  are  applicable  to  a winter  nighttime  propagation  model,  which  is  probably  the  most 
favorable  condition  for  the  operation  of  a communication  circuit.  We  then  consider  the 
increased  requirement  presented  by  the  summer  nighttime  noise,  and  finally  the  further 
increased  requirement  caused  by  lossy  cay  time  transioaospheric  propagation  during  both 
the  winter  and  summer. 
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Fig.  1— Illustration  of  propagation  path  from  ELF / 
VLF  transmitting  satellite  to  the  ionospheric  exit 
point.  The  Cartesian  coordinate  frame  for  the 
calculation  of  field  components  is  shown 


THEORY 

The  excitation  of  waveguide  modes  by  a satellite-borne  tiansmitter  has  been  consid- 
ered previously.  Galejs  [1,2]  gave  the  first  treatment  of  the  efficiency  by  which  an 
antenna  in  the  ionosphere  can  excite  a mode  in  ihe  earth-ionosphere  waveguide.  He 
showed  that  the  relative  effi  nency  of  a transmitting  antenna  in  the  ionosphere  is  at  most 
approximately  equal  to  that  of  an  optimally  located  antenna  in  the  earth  ionosphere 
waveguide,  and  it  is  usually  nuch  less  efficient.  His  numerical  results  are  for  frequencies 
between  10  and  1000  Hz  an'*  between  10  kHz  and  30  kHz  daytime  and  nighttime  propa- 
gation conditions.  Einaudi  and  Wait  [3,4]  give  some  numerical  results  for  a normalized 
excitation  factor  of  the  zero  o der  waveguide  mode  in  the  frequency  range  between  500 
Hz  and  7.0  kHz  for  the  case  of  a vertical  antenna  and  a vertical  magnetic  field.  The 
values  o ' the  relative  excitati  an  factor  are  much  smaller  than  those  predicted  for  an 
equivalent  transmitting  dipole  moment  optimally  oriented  within  the  tarth-ionosphere 
waveguide.  However,  Einaudi  and  Wait  assume  free  space  propagation  both  below  and 
above  the  plasma  slab  representing  the  D-region. 

Pappert  [5]  has  calculated  the  excitation  of  the  earth-ionosphere  waveguide  at  ex- 
vremely  low  frequencies  for  a satellite  at  9500-km  height  using  a realistic  nighttime  iono- 
sphere profile  and  several  geomagnetic  dip  angles  und  propagation  directions.  His  results 
indicate  that  at  75  Hz,  a horizontal  electric  d:pole  in  the  ionosphere  would  b?  approxi- 
mately 5 dB  more  efficient  than  a similar  horizontal  dipole  laid  on  a poorly  conducting 
(o  = 10'4 mhe/iV.)  earth  (similar  to  the  proposed  SANGUINE/SEAFARER  antenna).  How- 
ever, Pappert  states  that  an  “excessively  large”  magnetic  dipole  moment  would  be  re- 
quired to  give  long-range  communications  coverage  similar  to  a SANGUINE/SEAFARER 
transmitting  system. 
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A waveguide-mode  power  budget  at  a given  frequency  requires  that  the  mode  propa- 
gation parameters  be  known.  Barr  [7-11],  3alejs  [12],  and  Wait  [13]  have  given  numeri- 
cal values  of  the  waveguide-mode  propagation  parameters  in  the  frequency  range  of 
interest.  Barr  has  given  especially  complete  results  >n  the  1-10  kHz  range  for  numerous 
daytime  ionosphere  electron  density  profiles.  Computer  programs  [14]  are  available  for 
calculating  waveguide-mode  propagation  oarameters  for  arbitrary  electron  and  ion  density 
profiles,  geomagnetic  dip  angles,  and  propagation  direction". 

The  computer  program  of  Pappert,  etal.  [14]  determines  the  waveguide-mode  param- 
eters by  finding  the  roots  of  an  equation  involving  the  ionosphere  end  ground  reflection 
coefficients.  It  does  not  give  values  of  the  fields  within  the  ionosphv.."',  so  an  alternative 
approach,  such  as  that  of  Pitteway  [15] , must  be  used  to  determine  the  structure  of  the 
waveguide-mode  fields  within  the  ionosphere. 

In  the  present  study,  waveguide-mode  calculations  are  performed  in  the  frequency 
range  from  250  Hz  to  9 kHz.  Calculations  are  made  of  the  required  magnetic  moment 
which  a satellite  must  produce  to  provide  a usable  S/N  ratio  at  one  megameter  (Mm) 
from  the  ionospheric  exit  point.  The  power  needed  to  generate  this  magnetic  moment 
using  a loop  antenna  is  calculated  using  the  methods  of  Baker,  et  al.  [6] . Although  some 
effort  is  m?de  here  to  optimize  the  selection  of  satellite  height,  frequency  of  operation, 
and  loop  antenna  radius,  it  is  not  possible  to  optimize  the  system  fully  due  to  funding 
constraints.  A selection  of  alternate  frequencies  or  a different  loop  antenna  radius  might 
reduce  the  calculated  powers  significantly. 

Figure  1 shows  a Cartesian  cooidinate  system  in  which  the  vertical  coordinate  of  the 
transmitter’s  location  is  denoted  by  zs.  The  distance  between  the  point  on  the  earth’s 
surface  beneath  the  ionospheric  exit  point  and  the  receiver  is  y.  The  x axis  in  this  coordi- 
nate system  is  perpendicular  to  the  direction  of  wave  propagation.  An  arbitrary  dipping 
static  magnetic  field  can  be  inserted  in  this  system,  such  that  it  will  maintain  a constant 
angle  with  respect  to  the  direction  of  propagation  of  the  wave.  The  resultant  Maxwell’s 
equations  can  be  factored  into  z and  y-dependent  terms  and  solved  by  a separation  of 
variables  technique. 

The  relative  efficiency  rjjj'  \zg)  of  a magnetic  dipole  transmitter,  located  at  some 
height  zs  in  the  ionosphere,  for  exciting  the  nth  mode  of  the  earth-ionosphere  waveguide 
may  be  defined  as  the  ratio  of  a waveguide-mode  field  produced  by  a satellite  transmitting 
antenna  to  the  waveguide-mode  field  produced  by  an  equivalent  transmitting  antenna 
situated  on  the  earth’s  surface.  In  this  definition,  the  transmitting  antenna  on  the  earth’s 
surface  must  have  the  same  magnetic  dipole  moment  as  the  satellite-borne  antenna,  and 
be  located  at  a point  on  the  earth’s  surface  exactly  below  the  ionospheric  exit  point  of 
the  waves  transmitted  by  the  satellite.  The  waveguide-mode  fields  generated  by  the  two 
antennas  must  be  produced  at  the  earth’s  surface  at  the  same  distance  and  direction  from 
the  ionospheric  exit  point.  The  orientations  of  the  satellite  dipole  and  the  dipole  on  the 
earth’s  surface  must  be  specified  to  complete  this  definition.  Thus,  the  relative  efficiency 
Vj?'(zt)  of  a vertical  magnetic  dipole  at  height  zs  in  the  ionosphere  for  exciting  a given 
waveguide  mode  may  be  defined  as  the  ratio  of  the  z -component  of  the  electric  field  in 
'he  nth  waveguide  mode  produced  at  the  earth’s  surface  by  this  satellite  transmitting 
,m  tenna  to  tne  same  field  component  produced  by  an  equivalent  antenna  on  the  earth’s 
s irface.  In  this  defir.  ion,  the  transmitting  antenna  on  the  earth’s  surface  must  have  the 
srme  vertical  magnetic  dipole  moment  as  the  satellite-borne  antenna,  and  be  located  at  a 
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point  on  the  earth’s  surface  exactly  below  the  ionospheric  exit  point  of  the  waves  trans- 
mitted by  the  satellite.  The  waveguide-mode  folds  in  the  nth  waveguide  mode  must  be 
produced  at  the  earth’s  surface  at  the  same  distance  and  direction  from  the  ionospheric 
exit  point.  The  relative  efficiency  of  a hor.zontal  magnetic  dipole  transmitter  n the  iono- 
sphere can  be  defined  in  the  same  way;  however,  the  reference  ground-based  transmitting 
dipole  could  be  either  horizontal  or  vertical,  and,  if  horizontal,  it  could  be  oriented  either 
parallel,  perpendicular,  or  at  some  arbitrary  angle  to  the  direction  of  propagation  of  the 
waveguide  mode.  Because  of  the  flexibility  of  the  choice  of  the  reference  dipole’s  orienta- 
tion, there  are  at  least  nine  possible  efficiency  definition  combinations  for  magnetic 
dipole  transmitters,  and  nine  more  for  electric  dipole  transmitters.  We  will  restrict  our 
attention  in  this  paper  to  magnetic  dipole  antennas,  and  will  select  a reference  antenna 
with  the  best  orientation  for  exciting  the  mode  in  question. 

Galejs  [1]  and  Pappert  [5]  show  that  the  efficiency  ratio  as  discussed  above  is  equal 
to  the  ratio  of  the  respective  field  component  values  of  the  waveguide  mode  propagating 
in  the  opposite  direction  (from  the  receiver  toward  the  source),  where  the  numerator  of 
the  ratio  is  the  field  component  evaluated  in  the  ionosphere,  the  cenominator  is  the  field 
component  evaluated  at  the  earth’s  surface,  and  the  direction  of  the  earth’s  magnetic 
field  is  inverted.  Thus,  the  relative  efficiency  of  a vertical  electric  dipole  in  the  ionosphere 
to  one  on  the  earth’s  surface  for  exciting  a given  waveguide  mode  is  numerically  equal  to 
the  ratio  of  the  vertical  electric  field  present  in  that  waveguide  mode  at  that  height  in  the 
ionosphere  to  the  vertical  electric  field  carried  by  that  waveguide  mode  at  a point  on  the 
earth’s  surface  directly  below  the  ionospheric  exit  point  (when  the  direction  of  propaga- 
tion of  the  mode  and  the  direction  of  the  earth’s  magnetic  field  have  been  inverted). 

The  vertical  electric  field  Ez(0,  y)  generated  at  the  surface  of  the  earth  by  a trans- 
mitting satellite  loop  antenna  at  a height  zs,  having  current  / and  area  da,  can  be  calcu- 
lated from  the  following: 

E2(0,y)[by  Ida(zs)  l;)  Hj{zs)Ida(zs)  \j 

£,(0,y)[bytta(0)l,]  H;(0)lda(0)~  (1) 


where 

Ez(0,  y) [by  Ida(zs)  l;]  is  the  vertical  electric  field  at  the  earth’s  surface  gen  ’•ated  by 
a source  having  magnetic  dipole  moment  Ida  (zs)  at  the  height  z3,  oriented  with  the 
dipole  moment  along  the  y'th  axis, 

Ez( 0,  y)[by  /do(0)l(]  is  the  vertical  electric  field  r .he  earth’s  surface  generated  by 
a source  having  magnetic  dipole  moment  Ida( 0)  at  tne  earth’s  surface,  oriented  with 
the  dipole  moment  along  the  .th  axis, 

Hj(zs)  is  the  /th  component  of  the  magnetic  field  at  height  zs, 

Hj( 0)  is  ..he  ith  component  of  the  magnetic  field  at  the  surface  of  the  earth, 


A 


NRL  REPORT  8032 

lda(zs)\:  is  the  magnetic  dipole  moment  of  a loop  transmitter  located  at  height  ?s , 
orientea  with  the  dipole  moment  along  the  ;th  axis, 

I is  the  current  in  the  loop, 

da(zs)  is  the  area  of  the  loo^ 

Setting 


lda(zt)\,  « 7da(0)l,  (2) 

we  obtain 

E2(0,y)[by  Ida(zs)\j]  Hj(zs) 

r?v(2s)  " -E2tO,y)tby/aa(0)i,]  = ~H{0)  ’ (3) 

The  coordinates  i and  j can  be  arbitrarily  chosen  in  (1-3).  In  this  report  we  will  con- 
centrate on  the  TM-mode  fields,  because  they  are  most  effective  for  communications  to 
receivers  located  at  or  belov  the  earth’s  surface.  A magnetic  loop  transmitter  will  be  con- 
sidered because  the  power  losses  of  such  a transmitter  are  read;,y  calculable  [6].  The 
reference  magnetic  loop  transmitting  antenna  is  located  on  the  earth’1  surface  with  its 
axis  parallel  to  the  ,r  axis  (Fig.  1).  This  is  the  optimum  orientation  lor  generating  a TM- 
mode  wave  propagating  in  the  y direction,  and  is  most  meaningful  in  referring  the  TM- 
mode  fields  from  a satellite  transmitter  to  the  fields  produced  by  a loop  antenna  at  the 
earth’s  surface.  The  efficiency  of  a satellite  transmitter  is  then  given  by 


Hj{zs) 

Vxj(Zs)  " ~H^(0)  ‘ (4) 

Equations  for  the  fields  at  a point  at  distance  y from  a loop  transmitter  in  the  earth- 
ionosphere  waveguide  have  been  given  by  Oalejs  [12)  for  all  components  of  the  E and  H 
fields.  The  equations  for  the  contributions  by  the  TM-mode  fields  to  the  vertical  electric 
and  horizontal  magnetic  fields  are 

1<:z  (*,  ,v ) = co  r) 0I da  Y S™F„  ( 5 ) 
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where 


1 / y la 

F„  = — — J— G„(z)G„(z)A„exp(/7r/4  + ikaySn)  (7) 

V si.i  (y/a) 

a = the  radius  of  the  earth  (m), 
z = the  receiver  height  (m), 
zs  = the  transmitter  height  (m), 

A„  = the  excitation  factor  of  the  nth  mode, 

S„  = the  sine  of  the  eigenangle  of  the  nth  waveguide  mode, 
y = the  great  circle  distance  between  the  wave  exit  point  and  receiver  (m), 

X = the  free-space  wavelength  of  the  wave  (m), 
k0  - 27T/X  (m*1), 

1 = current  of  loop  (amps), 

G„  = the  height  gain  function  of  the  nth  mode, 

da  = area  of  loop  (m2', 

co  = angular  frequency  of  wave  = 2 vf  Hz, 

H0  = permeability  of  free  space  = 4tt  X 10'7, 
h = the  height  of  the  waveguide  (in). 

Equations  (5-7)  apply  to  the  case  in  „hich  both  the  transmitter  and  receiver  are 
within  the  waveguide  (0  < z,  zs  < h)  The  sine  of  the  mode  eigenangle  S„  is  related  to 
the  attenuation  rate  a and  the  phase  velocity  vph  of  the  waceguide-mode  waves  accord- 
ing to  Galejs  (12| : 


Re  (S„)  = CjVp 

(8) 

a = 0.02895  co!m(Sn)(in  dli/Mm). 

(9) 

In  (8),  c is  the  speed  of  light. 

The  following  procedure  is  used  to  calculate  i?,;(cs).  First,  the  mode  eigenangle  and 
the  vai "“s  of  the  electric  fie!.!  components  at  several  heights  in  the  waveguide  are  calcu- 
lated according  to  the  mi  mod  of  I’appert,  et  al.  jl4|.  Based  on  the  eigenangle  of  an  inci 
dent  wave  on  the  lower  ionosphere,  the  wovclield  of  the  penetrating  and  nonpenctrating 
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waves  are  calculated  at  the  same  heights  in  the  waveguide  and  at  selected  heights  in  the 
ionosphere  according  to  the  procedure  of  [15-17], 

Next,  the  coefficients  bn  and  up  are  determined  applying  Cramer’s  Rule  to  (10-12). 
In  this  way,  the  waveguide  mode  may  be  constructed  from  a linear  combination  of  the 
penetrating  and  nonpenetrating  field  solutions: 


Ez^n\z)  = bpE%en(z)  + bnEzon{z) 

E<yn\z)  = bpE?n(z)  + LnEnyon(z) 
Elxn\z)  = hpE*”(z)  + bnE™'\z). 


In  the  above  equations, 

E^n\z),  E^n\z),  and  E^zn\z)  are  the  components  of  the  nth  waveguide-mode  elec- 
tric vector  at  height  z in  the  waveguide, 

E%en(z),  E?,en(z),  and  E,’en(z)  are  tne  components  of  the  electric  vector  of  the 
penetrating  wave  at  the  height  z, 

E"on(z),  Eyon(z),  and  E'Jon(z)  r. re  the  components  >f  the  electric  vector  of  the  non- 
penetrating wave  at  the  height  z. 

Coefficients  bp  and  bn  giving  the  contributions  of  the  penetrating  and  nonpenetrat- 
ing waves  in  forming  a waveguide  mode.  Any  pair  of  (10-12)  constitute  two  equat.o.is  in 
two  unknowns  (bp,  btl),  and  permit  the  values  of  bp  and  bn  to  be  readily  determined  by 
applying  Cramer’s  Rule.  The  consistency  of  the  values  of  the  third  equatnn  may  be  used 
as  an  independent  check.  In  addition,  similar  sets  of  field  values  may  be  obtained  at  other 
heights  in  the  waveguide  and  additional  sets  of  equations  used  to  check  the  accuracy  of 
bn  and  bp.  The  wave  gu'de-mode  fields  in  the  ionosphere  may  be  readily  determined  from 
the  and  bp  amplitudes  and  the  values  of  the  penetrating  and  nonpenetrating  wave- 
field:.  in  the  ionosphere.  The  efficiency  n[!y(zs)  for  the  nth-order  TM  mode  is  given  by 


4nj](zt) 


yr  + bniijion(zs 
bpHPcn(  0)  + b, ,//;!""( 0) 


where  ll^en(zs)  and  //"0,,(2S)  are  the  magnetic  fields  of  the  penetrating  and  nonpenetratir.g 
waves  at  height  zs. 

To  determine  the  waveguide-mode  fields  at  some  distance  y from  the  ionospheric 
exit  point,  we  multiplied  the  terms  m the  mode  summations  of  the  right-hand  side  of 
(5-6)  by  the  efficiency  value  ijxyUs ) of  the  appropriate  mode  for  the  appropriate  satellite 
loop  orientation  and  height.  The  transmitter  height  gain  function  C„(z..)  which  enters  in 
the  calculation  of  En  is  replaced  by  the  function  n['j\zs)-  The  results  are 
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f ■ 1 ^ )/M^)G"iz)4"iH*i}An  exp(“/4 + ik°ys) ■ (i4) 

The  approximation  of  geometric  optics  may  be  applied  to  obtain  an  approximation 
to  the  wave  guide-mode  fields  at  considerable  heights  above  the  lower  D and  ^-ionospheric 
region.  The  equations  of  Ginzburg  [18]  may  be  expressed  as 


Ex,y(z)  = {Exy)0(nQ\n(z))xl2exp 
Hx,y(z)  = (^,y)o(n(z)i«o)1/2  exP 


* 

> 

l 

o 

t 

fz  V 

iii; 

| n(z)dz\ 
z0  1. 

r : 

t - cx  | 

~2  V 

iiv 

n(z)dz  \ 

\ 

J 

\ 

!0  I. 

(15) 

(16) 


(Exy)0  and  (Hx  y)0  are  electric  and  magnetic  fields  at  height  z0, 

Ex  y(z)  and  Hx  y(z)  are  electric  and  magnetic  fields  at  height  z, 

n0  and  n(z)  are  the  indices  of  refraction  of  the  propagating  wave  at  heights  z0  and 
z,  respectively. 

Thus,  the  magnetic  field  achieves  a maximum  value  where  the  index  of  refraction  is 
maximum.  The  value  of  the  tj(.")(zs)  efficiency  ratio  may  be  written 


4-'/' ) (2s ) vinj  *(  1 20/em ) [r?  (z  ) |rj(  120fe  m ) ] 1 12 


IXJ 

r z 

[X  exp[i'u»((  - c'1  ri(z)dz] 


(17) 


in  terms  of  the  efficiency  at  120  km,  the  indices  of  refraction  at  z0  (120  km),  and  the 
transmitter  height  zs. 

The  input  power  to  the  antenna  required  to  produce  a 0-dB  S/N  ratio  in  a 1-Hz 
bandwidth  is  given  by 


^req  ^req  ( R 0|,  m + R ra(j ) 


(18) 


where 


/?ohm  *s  ohmic  resistance  of  the  loop, 
Rtacj  is  the  radiation  resistance  of  the  loop, 


(19) 
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N(f)  = the  rms  effective  magnetic  field  noise  in  a 1-Hz  bandwidth, 
da  - the  area  of  the  loop  antenna. 

In  the  above  equation,  t)\"\zs)  is  given  by  (13)  between  68  and  120  km,  and  by  (17) 
for  transmitter  heights  above  120  km. 


CALCULATIONS  AND  RESULTS 

For  nighttime  propagation  conditions  the  Deeks  wr^r  nighttime  Ionosphere  electron 
density  profile,  as  given  by  Galejs  (12],  was  used  (Fig.  2)  with  the  value  of  collision 
frequency  as  given  below: 


10"  10°  :0'  IC?  I05  I04 

ELECTRON  DENbiTY  Ne-cm> 

Fig.  2— Ionospheric  electron  density  profile  (From 
Fig.  22,  Terrestrial  Propagation  of  Long  Elec- 
tromagnetic Waves,  by  J.  Galejs,  Peragmon 
Press,  Copyright  1972  by  J.  Galejs;  used  by  per- 
mission.) 


vJz)  = 1.8158  X 1011  exp  (-  15z). 


Figures  3,  4,  and  5 show  propagation  parameters  as  a function  of  frequency  for  the 
Deeks  winter  nighttime  electron  density  profile.  Results  are  given  for  three  different  wave- 
guide modes  that  propagate  readily  between  0.5  and  10.0  kHz.  The  direction  of  propaga- 
tion for  this  mode  was  due  north  and  the  dip  angle  was  60°. 

Figures  6 and  7 show  the  height  variation  of  the  efficiency  ratio  for  Quasi-TE  and 
Quasi-TM  mode  fields  at  6 kHz.  Figure  8 shows  the  efficiency  variation  for  the  zero  order 
mode  at  250  Hz. 


Figures  9 and  10  give  the  efficiency  ratio  for  9-kllz  Quasi-TM  and  Quasi-TE  modes. 
In  calculating  the  results  of  Figs.  3-10,  the  ionospheric  electron  density  parameters 
between  100  and  120  km  were  held  constant,  equal  to  the  density  at  100  km. 
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Fig.  3— Lowest-ordcr  TEM  mode  propagation  parameters  vs  frequency 


To  determine  r^x'j\zs)  for  heights  between  100  and  1000  km,  a ray  trace  was  per- 
formed using  the  procedures  and  ionosphere  model  A-3  (for  a winter  nighttime  iono- 
sphere) described  by  Kelly,  et  al.  [19].  Figure  11  shows  curves  of  «(;:)  for  frequencies 
between  250  Hz  and  10  kHz. 

The  net  radiation  resistance  of  a loop  antenna  having  a radius  of  100  m was  calcu- 
lated according  to  procedures  given  by  Baker,  et  al.  16]  as  a function  of  ionospheric 
height  at  frequencies  between  250  Hz  and  10  kHz,  and  is  displayed  in  Fig.  Vl.  Results 
are  given  for  tve  loop  oriented  parallel  or  perpendicular  to  the  magnet:"  field.  Ohmic 
resistances  of  loop  antennas  are  readily  calculated.  For  cases  in  which  the  ohmic  loss  is 
much  greater  than  the  radiated  power,  it  is  usually  possible  to  improve  the  ratio  of 
radiated  power  to  ohmic  power  by  increasing  the  loop  radius.  At  250  Hz  and  1 kHz.  we 
used  the  zero-c  rder  TEM-mode  parameters.  At  6000  and  9000  kHz,  we  used  the  modal 
parameters  of  the  first-order  TM  mode.  Additional  modes  would  certainly  contribute  to 
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the  fields  at  the  earth's  surface  at  1900  km  from  the  subsatellite  point,  but  the  dominant 
contribution  would  normally  be  provided  by  these  two  low-order  modes. 


Tables  1 througi  \ present  the  required  satellite  transmitter  powers  for  severrl  trans- 
mitter frequencies,  and  altitudes  for  various  noise  ar.d  propagation  conditions.  In  these 
tables  the  radiated  power,  ohmic  loss  power,  and  total  power  are  denoted  by  Pra(l, 

P ohmic ' anc*  ? total’  respectively.  Table  1 gives  ‘N  results  calculated  for  winter  nighttime 
noise  and  propagation  condition.!.  Tables  2,  3,  and  4 contain  the  results  for  summer  day 
and  night  and  for  winter  day  conditions,  and  are  simply  an  extrapolation  of  the  winter 
nighttime  results  by  the  application  of  appropriate  propagation  loss  and  atmospheric 
noise  differences.  This  simplification  of  the  calculations  is  justified  in  light  of  the  un 
certainties  and  variability  of  the  ionospheric  models.  Thus,  Table  2 was  generated  for 
summer  nighttime  conditions  by  using  the  same  propagation  calculations  as  in  Table  1, 
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Fig.  5— First-order  TE  mode  propagation  parameters  vs  frequency 


and  by  increasing  the  atmospheric  noise  level  as  predicted  by  Watt  [20]  for  the  summer 
nighttime  Mediterranean  Sea  area.  The  values  of  the  daytime  and  nighttime  atmospheric 
noise  at  250  Hz  in  the  Mediterranean  Sea  are  taken  as  -140  dB  below  1 amp/m  ’ 'ring 
summer  and  -146  dB  below  1 amp/m  during  the  winter  period.  These  values  are  charac- 
teristic of  75-Hz  measured  v-lues  [21]  and  are  probably  higher  than  expected  at  250  Hz. 

At  6 kHz  the  efficiency  factor  r)xx  (120)  is  approximately  6 dB  lower  for  a wave 
propagating  south  from  the  ionospheric  exit  point  than  for  a wave  propagating  north.  At 
250  Hz  and  1 kHz  the  efficiencies  for  exciting  the  north  and  south  propagating  waves 
differ  by  less  than  2 dB.  The  calculated  results  in  Tables  1 through  4 are  valid  for  tire 
case  in  which  the  wave  propagates  north  from  the  ionospheric  exit  point. 

The  predictions  of  Tables  3 and  4 for  summer  and  winter  daytime  conditions  were 
obtained  using  Watt’s  [20]  noise  values  for  the  appropriate  times  and  by  increasing  the 
transionospheric  propagation  loss  by  the  amounts  derived  from  Table  5.  The  nighttime 
attenuations  sho'  ->  in  Table  5 were  deduced  by  averaging  separately  the  daytime  and 
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nighttime  loss  values  predicted  by  the  calculations  of  Pitteway  and  Jespersen  [15], 
Booker,  Crain,  and  Field  [22],  Al’pert  and  Fligel  [23],  and  Helliwell  [24].  The  day- 
time power  requirements  were  also  increased  by  an  additional  factor  because  of  the  extra 
attenuation  that  the  waves  undergo  in  propagating  from  the  ionospheric  exit  point  to  the 
receiver  location  at  1 Mm  range.  The  calculations  of  Barr  [7]  for  six  different  ionospheric 
modeis  were  averaged  to  determine  « mean  daytime  attenuation  rate  for  the  TFM  and  TM 
modes  used  in  the  present  calculatio,  is.  The  average  attenuation  rates  are  given  in  Table  6. 


Fig.  8— Relative  excitation  efficiencies  vs  height  of  the  zero-order  TEM 
mode  at  0.25  kHz 
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9 kHz  OUASI-TE  MODE  FIELDS 


Fig.  9 Relative  excitation  efficiency  vs  height  of  the  first 
order  TK  mode  at  9.0  kHz 
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RADIATION  RESISTANCE  (OHMS) 

Fig.  12— Radiation  resistances  of  a 100-m  radius  loop  at 
several  frequencies  and  heights  in  the  ionosphere,  and  for 
loop  dipole  orientations  parallel  (Rj)  and  perpendicular 
(/?l)  to  the  earth’s  magnetic  field.  Values  of  R j were  used 
jf  j in  computing  Tables  1 through  4. 
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Table  1 — Compilation  of  Required  Power  Results  for  a 104  -kg,  100-m 
Radius  Loop  Transmitting  Antenna  for  y = 1 Mm* 


zs(km) 

^rad(kW) 

^ohrnic(kW) 

^total(kW) 

120 

1.3  X 101 

1.9  X 104 

1.9  X 104 

300 

1.2  X 102 

3.6  X 103 

3.7  X 103 

800 

1.8  X 104 

9.4  X 103 

2.7  X 104 

120 

4.0  X 10° 

1.3  X 102 

1.4  X 102 

300 

7.4  X 10; 

2.7  X 101 

1.0  X 102 

800 

2.8  X 10° 

7.3  X 101 

7.6  X 101 

120 

1.6  X 10° 

8.6  X 10“2 

1.7  X 10° 

300 

1.8  X 101 

1.5  X 10~2 

1.7  X 101 

800 

6.4  X 10'1 

4.0  X 10'2 

6.7  X 10'1 

120 

4.5  X 101 

1.1  X 10° 

4.7  X 10 1 

300 

1.5  X 101 

2.5  X 10'1 

1.5  X 101 

800 

4.5  X 101 

6.8  X 10’1 

4.5  X 101 

♦Winter  nighttime  propagation  and  noise  with  S/N  = 0 dB  in  a 1-Hz  bandwidth  (Mediterranean  Sea  area). 


Table  2 — Compilation  of  Required  Power  Results  for  a 104  -kg,  100-m 
Radius  Loop  Transmitting  Antenna  for  y = 1 Mm* 


A Hz) 

zs(km) 

^rad(kW) 

^ohmie(kW) 

250 

120 

5.2  X 101 

7.7  X 104 

300 

4.7  X 102 

1.4  X 104 

800 

7.2  X 104 

3.7  X 104 

1000 

120 

5.0  X 101 

1.7  X 103 

300 

9.3  X 102 

3.4  X 102 

800 

3.5  X 101 

9.2  X 102 

6000 

120 

2.6  X 101 

1.1  X 10° 

300 

2.8  X 102 

2.4  X 10*1 

800 

1.0  X 101 

6.6  X 10*1 

9000 

120 

1.8  X 102 

4.5  X 10° 

300 

6.0  X 102 

1.0  X 101 

800 

1 8 X 102 

2.7  X 10° 

JWW) 

7.7  X 104 

1.5  X 104 
1.1  X 105 

1.7  X 103 
1.3  X 103 

9.5  X 102 

2.7  X 101 

2.8  X 102 
1.1  X 101 

1.8  X 102 
6.0  X 102 

1.8  X 102 


♦Winter  nighttime  propagation  and  summer  nighttime  noise  with  S/N  = 0 dll  in  a I Hz  bandwidth 
(Mediterranean  Sea  area). 
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Table  3 — Compilation  of  Required  Power  Results  for  a 104  -kg,  100-m 
Radius  Loop  Transmitting  Antenna  for  y = 1 Mm* 


A Hz)  zs(km) 

^rad(kW) 

^ohmic(kW) 

Ptotal(kW) 

8.2  X 10J 

7.4  X 103 

1.1  X 108 

5.4  X 104 

1.0  X 108 

3.8  X 104 

3.1  X 104 

3.3  X 105 

1.2  X 104 

1.1  X 105 

3.8  X 10 5 

1.1  X 105 


1.2  X 106 

2.3  X 108 

5.9  X 105 

1.8  X 106 

3.7  X 105 

9.9  X 105 

1.3  X 103 

2.8  X 102 

7.8  X 102 

3.0  X 103 

6.5  X 102 

1.7  X 103 


1.2  X 108 

2.4  X 108 

1.7  X 108 

1.9  X 106 

1.4  X 108 

1.0  X 108 

3.2  X 104 

3.3  X 105 

1.3  X 104 

1.1  X 105 

3.8  X 105 

1.1  X 105 


Daytime  propagation  and  summer  daytime  with  S/N  = 0 dB  in  a l*Hz  bandwidth  (Mediterranean  Sea  area). 


Table  4 — Compilation  of  Required  Power  Results  for  a 104-kg,  100-m 
Radius  Loop  Transmitting  Antenna  for  y = 1 Mm* 


z,  (km) 


^rad(kW) 

^ohmic(kW) 

WkW) 

2.1  X 102 

1.9  X 103 

2.9  X 105 

3.1  X 105 
5.7  X 104 
1.5  X 105 

3.1  X 105 
5.9  X 104 
4.3  X 105 

8.3  X 102 
1.5  X 104 
5.8  X 102 

2.8  X 104 
5.7  X 103 
1.5  X 104 

2.9  X 104 
2.1  X 104 
1.'  X 104 

5 8 X 10- 

6.2  X 103 

2.3  X 102 

2.6  X 101 

5.4  X 10° 

1.5  X 101 

6.0  X 102 
6.2  X 103 
2.4  X 102 

9.8  X 103 
o.3  X 104 

9.8  X 103 

2.5  X 102 

5.5  X 101 
1.4  X 101 

9.8  X 103 
3.3  X 104 
9.8  X 103 

Daytime  propagation  and  winter  daytime  noise  with  S/N  = 0 dB  in  a 1-Hz  bandwidth  (Mediterranean  Sea 
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Table  5 — Nominal  Ionospheric  Transmission  Losses  Derived  From 
the  Calculations  of  [15]  and  [22-24] 


Avg.  Daytime 

Std.  Deviation 

Loss  (dB) 

Daytime  Models  (dEi) 

9.8 

3.2 

9.2 

2.1 

16.8 

3.9 

21.4 

5.7 

Avg.  Nighttime 
Loss  (dB) 


Std.  Deviation 
Nighttime  Models  (dB) 


Table  6 - Nominal  Daytime  Attenuation  Rates  in  the  Sarth-Ionosphere 
Cavity  Based  on  Six  Models  Used  by  Barr  [7] 


f(  Hz) 

Avg.  Daytime 
Attenuated  Rate  (dB/Mm) 

Std.  Deviation 
Daytime  Models  (dB/Mm) 

250 

3.8 

1000 

16.2 

6.1 

6000 

8.5 

3.3 

9000 

4.3 

0.5 

DISCUSSION 


The  cases  analyzed  in  Tables  1-4  should  not  be  considered  as  selected  optimum  con- 
figurations or  frequencies.  The  calculations  only  illustrate  the  broad  range  of  required 
transmitter  powers  that  can  be  obtained  for  different  communication  link  requirements. 

For  the  cases  studied  here,  the  time  during  which  the  satellite  and  groundbased  re- 
ceiver are  in  contact  with  each  other  is  relatively  brief.  Assuming  that  the  ionospheric 
exit  point  is  moving  at  a typical  satellite  velocity  of  5 km/s,  the  maximum  contact  time 
vhiie  a 1000-km  radius  coverage  area  moves  over  a fixed  receiver  location  would  be  ap- 
proximately 7 min. 

In  the  results  shown  in  Tables  1-4,  we  have  not  optimizcu  the  transmitting  loop 
antenna  radius  to  the  full  extent  possible.  Further  calculations  indicate  that  at  the  250-llz 
frequency  considered  here,  a more  balanced  distribution  of  input  power  between  radia- 
tion and  ohmic  losses  may  be  accomplished  by  increasing  the  loop  antenna  radius  beyond 
the  value  of  100  m assumed  for  ".'able  1.  This  increase  of  radius  (keeping  the  total  mass 
fixed)  tends  to  increase  the  ohmic  resistance  as  the  square  of  the  loop  radius.  The  loop 
radiation  resistance  (and  the  power  injected  into  the  earth-ionosphere  waveguide)  in- 
creases in  proportion  to  ihe  fourth  power  of  the  loop  radius.  For  example,  using  a 
1000-rn  radius  loop,  this  type  of  extrapolation  implies  that  under  winter  nighttime  con- 
ditions, the  total  required  power  for  a 250-Hz  satellite  transmitter  at  a 300-km  altitude 
height  is  614  kW  with  470  kW  radiated  and  144  kW  of  ohmic  losses  using  a 10‘'-kg 
antenna.  A more  exact  calculation  at  250  Hz  leads  to  the  total  power  requirement  of 
799  kW  with  736  kW  radiated  and  63  kW  of  ohmic  losses.  Obviously,  some  provision 
must  be  made  for  radiating  the  heat  generated  by  the  ohmk  losses  of  the  antenna. 
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The  results  also  indicate  that  the  transmitter  power  required  to  produce  a given  S/N 
ratio  at  a given  point  on  the  earth’s  surface  ranges  from  less  than  a kilowatt  to  thousands 
of  megawatts,  depending  on  the  many  variables  affecting  the  signal  and  noise  levels.  If  we 
limit  ourselves  to  ihe  results  for  a 6-kIlz  transmitter  operating  at  800  kni,  tne  power  re- 
quired to  produce  a 0-dB  S/N  ratio  in  a 1-Hz  bandwidth  varies  from  less  than  a kilowatt 
for  winter  lighttime  conditions  to  several  megawatts  for  summer  davtime  comb* ions.  ThL 
large  variation  occurs  because  the  noise  varies  by  about  19  dB  and  the  signal  by  about  13 
dB.  In  addition,  the  poorest  signal  and  highest  noise  occur  during  summer  day,  while  the 
best  signal  and  lowest  noise  occur  during  winter  night. 
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One  may  conclude  from  this  study  tliat  it  is  certainly  possible  to  achieve  some  com- 
munications during  the  optimum  period  using  relatively  low-powered  (few  kilowatts) 
transmissions,  but  that,  at  the  sample  frequencies  considered  (250  Hz,  1 kHz,  and  9 kHz), 
a truly  effective  system  would  require  increasing  powers  up  to  the  several-megawatt 
range  and  beyond,  since  greater  data  rate  and  coverage  area  are  always  desirable. 


The  U.S.  National  Aeronautics  and  Space  Administration  (NASA)  and  the  Energy 
Research  and  Development  Administration  (ERDA)  are  studying  plans  and  proposals  for 
future  space  nuclear  and  solar  electric  power  sources  which  range  in  power  from  a few 
kilowatts  to  thousands  ut  megawatts  [25-33].  The  expected  operation  of  the  reusable 
NASA  space  shuttle  would  make  the  launch  and  deployment  of  a reasonably  large  power 
source  and  transmitting  antenna  feasible  and  possibly  economical. 
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